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Abstract 

A Bianchi type-I cosmological model in the presence of a magnetic flux 
along a cosmological string is investigated. The objective of this study is to 
generate solutions to the Einstein equations using a few tractable assump- 
tions usually accepted in the literature. The analytical solutions are supple- 
mented with numerical and qualitative analysis. In the frame of the present 
model the evolution of the Universe and other physical aspects are discussed. 
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1 Introduction 



Since the observed Universe is almost homogeneous and isotropic, the space-time 
is usually described by a Friedman-Lamaitre-Robertson- Walker (FLRW) cosmol- 
ogy. Although this is a good approximation, the recent observations suggest the 
necessity of exploring beyond it. Also it is believed that in the early Universe the 
FLRW model does not give a correct matter description. The anomalies found 
in the cosmic microwave background (CMB) and the large structure observations 
stimulated a growing interest in anisotropic cosmological models of the Universe. 

Recently, cosmic strings have drawn considerable interest among the resear- 
chers for various aspects such as the study of the early Universe. The presence of 
cosmic strings in the early Universe could be explained using grand unified theo- 
ries. These strings arise during the phase transition after the Big Bang explosion 
as the temperature goes down below some critical point. It is believed that the ex- 
istence of strings in the early Universe gives rise to the density fluctuations which 
leads to the formation of the galaxies. Also the cosmic strings have been used in 
attempts to investigate anisotropic dark energy component including a coupling 
between dark energy and a perfect fluid (dark matter) [[T]|. The cosmic string is 
characterized by a stress energy tensor and it is coupled to the gravitational field. 

The presence of magnetic fields in galactic and intergalactic spaces is well 
established and their importance in astrophysics is generally acknowledged (see 
e. g. the reviews [2, 3] and references therein). In spite of the fact that the 
present day magnitude of the magnetic energy is very small in comparison with 
the matter density, it might not have been negligible during early stages of the 
evolution of the Universe. Any cosmological model which contains magnetic 
fields is necessarily anisotropic taking into account that the magnetic field vector 
implies a preferred spatial direction. 

Among many possible alternatives, the simplest and most theoretically ap- 
pealing of anisotropic models are Bianchi type-I (BI). For studying the effects of 
anisotropics in the early Universe on the present observed Universe, BI models 
have been investigated from different points of view. In this class of models it 
is possible to accommodate the presence of cosmic strings as an example of an 
anisotropy of space-times generated by one dimensional topological defects. 

In what follows we shall investigate the evolution of BI cosmological models 
in presence of a cosmic string and magnetic fluid. The paper has the following 
structure. We shall review the basic equations of an anisotropic BI model in the 
presence of a system of cosmic string and magnetic field. In Section 3 we intro- 
duce a few plausible assumptions and investigate their outcomes. The objective 
of this treatment is to generate exact solutions to the Einstein equations, supple- 
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mented with numerical and qualitative analysis. At the end we shall summarize 
the results and outline future prospects. 



2 Fundamental Equations and general solutions 

The line element of a BI Universe is 

ds^ = {dtf - ai {tfidx^f - a2{tf{dx^f - {tf{dx^f . (2.1) 

There are three scale factors a, (? = 1,2,3) which are functions of time t only and 
consequently three expansion rates. In principle all these scale factors could be 
different and it is useful to express the mean expansion rate in terms of the average 
Hubble rate: 

H=Uil + "^ + "^), (2.2) 

where over-dot means differentiation with respect to t. 

In the absence of a cosmological constant, the Einstein's gravitational field 
equation has the form 





= ktI, 


(2.3a) 


a2 CI3 0,2 CI3 




Ht. d] di, d\ 
—+—+—— 




(2.3b) 


^3 a\ ci\ 




d\ a? d\ d2 
—+—+—— 


= 


(2.3c) 


a\ a2 a\ a,2 






ai a2 1 <5[2 6(3 1 Q'i a\ 




(2.3d) 


a\ a2 0.2 C13 0.2 a\ 





where K is the gravitational constant. The energy momentum tensor for a system 
of cosmic string and magnetic field in a comoving coordinate is given by 

= pu^u" - Xx^x" +El , (2.4) 

where p is the rest energy density of strings with massive particles attached to 
them and can be expressed as p = + A , where pp is the rest energy density of 
the particles attached to the strings and A is the tension density of the system of 
strings [4", '5^, 165 which may be positive or negative. Here ut is the four velocity and 
Xi is the direction of the string, obeying the relations 

m/m' = —XiX^ = 1, Uix' = 0. (2.5) 
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In (12.41 ) E'^v is the electromagnetic field given by Lichnerowich [T|. In our 
case the electromagnetic field tensor F"^ has only one non- vanishing component, 
namely 

F23=/i, (2.6) 

where h is assumed to be constant. For the electromagnetic field Ej^ one gets the 
following non-trivial components 

1 2 3 

Eq^Ei ^ -E2 = -£3 = 2 2 • (^•'^^ 

2jUa2a3 

where /i is a constant characteristic of the medium and called the magnetic per- 
meability. Typically p. differs from unity only by a few parts in 10^ (/i > 1 for 
paramagnetic substances and p. < I for diamagnetic). 

Choosing the string along direction and using comoving coordinates we 
have the following components of energy momentum tensor [SJ: 

r- ^ A + (2.8b) 
Ti - -0. (2.80 
7^3^ = (2.8d) 

2/i 

where we introduce the volume scale of the BI space-time 

T = aiaja^ , (2.9) 

namely, T = [|9l. It is interesting to note that the evolution in time of T is 
connected with the Hubble rate (12.101) : 

- = 3H. (2.10) 

T 

In view of = from (llJbl) . (l23cl) one finds 

a2 = C3a3exp(^d J , (2.11) 
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with C3, J some real constants. 

Since the metric functions can be expressed in terms of T, let us first derive 
the equation for T. Summation of Einstein equations (I2.3al) . (|2.3bl) . (|2.3cl) and 3 
times (|2.3dl) gives 

T 3 / X a^' 



Taking into account the conservation of the energy-momentum tensor, i.e., 
T^.y = 0, after a little manipulation of (|2.8I) one obtains 

P + -P-— A=0. (2.13) 



3 Some examples and explicit solutions 

The above equations involve some unknowns and we need some supplementary 
relations between them to have a tractable problem. It is customary to assume a 
relation between p and A in accordance with the state equations for strings. The 
simplest one is a proportionality relation lH: 

p = aX. (3.14) 

The most usual choices of the constant a are 

{1 geometric string 

1 + CO CO > 0, p string or Takabayasi string (3.15) 
— 1 Reddy string. 

In order to solve the Einstein equations completely, we need also to impose 
some additional conditions. In what follows we shall illustrate the general consid- 
erations by two examples. 



3.1 Case 1: Hubble rate proportional to an eigenvalue of the 
shear tensor 

As a first example, we shall follow the condition introduced by Bali [TO] assuming 
that the average Hubble rate H (12.21) in the model is proportional to the eigenvalue 
a/ of the shear tensor o^. We shall only give the briefest account here and for 
further information the reader should consult [fTTIl . 
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For the BI space-time we have 



^l^-lU'ii+'^+'h.), (3.16) 



3 V ai a2 6(3 
Writing the aforementioned proportionality condition as 

H=pal, (3.17) 

one comes to the following relation 

fll = C23(«2«3)^, (3.18) 

where q= —{p+l)/ {4p + 1) is related to the proportionality constant p and C23 
is an integration constant. 

In this case (12.131) takes the form 

P+ip-^^)l = ^^ (3.19) 
^ q+l T 

and (12.121) now reads 

f = lK{p + |)t + Z7it('?-i)/(^+i) , (3.20) 

where bi = ^^c^^^'^^^^ is another constant. 

Let us now study (13.191) and (13.201) for different equations of state. Assuming 
the relation (13.141 ) between the pressure of the perfect fluid p and the tension 
density A from (13.191) one finds 

^ = (^^-1)-, (3.21) 

with the solution 



p=Z72T«(?+i) , (3.22) 

while the equation for t reads 

3 1 1 ?-i 

f = -k:Z?2(1 + — )t«(9+i) +Z?iT?+i , (3.23) 

^ 3 (JC 

involving another constant of integration b2. 

This equation can be set in the following form 

'Kb2{3a + l){g+l) ^ Mj+n ,2,/(,+i) + ^3 (3.24) 

q + a{q+l) q 

where b^, is an integration constant. 

More details for this model characterized by (|3.17l) including some numerical 
computations are given in flllll . 



6 



3.2 Case 2: Constraints on relative shear anisotropy parame- 
ters 

The second case is more involved taking into account that the anisotropy is con- 
nected with the values of the shear distortions. Besides the generalized Hubble 
parameter (12.21 ) we consider two relative shear anisotropy parameters defined by 
a2j|: 

R 

S 

When R = S = the universe will be the isotropic flat Friedman Universe. In this 
second example we shall assume a deviation from the Friedman model. 

We shall consider three different assumptions regarding the relative shear ani- 
sotropy parameters. Although these assumptions are distinct, they entail mild time 
dependencies for R and S. More exactly the time evolution of the scale factors 
is assumed to be very similar to that of the average Hubble rate H. 

3.2.1 Assumption 1: i? is constant 

Let us assume that R — ri with ri a constant. In view of (12.21) one finds, 

ai =C2a2T'''/^ (3.26) 

Then, together with (12.91 ). (12.1 II) we find the following expressions for metric 
functions 

ai = (c2c3)i/3r(3+2n)/9,f/f ^ (3.27a) 

a2 = (c2lc3)l/3T(3-'-i)/9,f/f ^ (3_27b) 

a, = (c2ci)-i/3T(3-n)/9,^/f . (3.27c) 
In this case (|2.13l) takes the form 

P+(p-^^)^-|^=0. (3-28) 
whereas, for t in this case we have 




(3.29) 
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) 



1 ,di _ «3\ 
H ai as 



(3.25a) 
(3.25b) 



Comparing this equation with the corresponding one (13.231) from Case 1, it 
resuks that in Case 2 the equation of evolution for T is more intricate. In ifTTIl 
we gave some asymptotic solutions of this equation. Here we give more general 
solution to the equation in question. 

In what follows we express equations (|3.28l) and (13.291) as a new system of 
differential equations which is easier to analyze: 

t = 3//T, (3.30a) 
f = -, (3.30b) 



T 



3aT 3a ' 

This system was investigated qualitatively using numerical methods. To nu- 
merically integrate the above system of differential equations with rational poly- 
nomials in the right hand sides we use the following method: each of the variables 
V = {x,H,p,T} was substituted by V = Vs/Vc with + = \. Further we find 
the common denominator of the right hand sides and moving away it we find the 
new system of equation with new parameters. This removing is equivalent to the 
substitution of an independent variable. In the new variable the right hand sides of 
the differential equations, as well as the integrable functions themselves, possess 
finite variation. So in the compact region the integral curves happen to be stable 
if the initial value is given inside the region. It does not change if we numerically 
integrate them in natural parameters (along the length of the curve). 

For exemplification we make the following specific choices: ri = 5.25, d ~ 
1.5, K = I, p. = I, 3 = 1,C3 = 1, and C2 = 1. We have investigated the system for 
two different values of a, namely, a = 1 and a = — 1 . To plot the graph^l it is 
more convenient to use instead of the function /, arctan/ in order to convert their 
infinite ranges of variation into finite ones. The initial condition we give on the 
line where H = and study the evolution of all the four functions both through the 
past and future. In the Figs. [T]and|2]we plot the evolution of T for a positive and 
negative a, respectively, while in Figs. [3]and|4]we do the same for energy density. 
In Figs. [5]and[6l 3D graph of T, H and p has been illustrated for a positive and 
negative a, respectively for the following initial values: 



Figures of this article are in color in the electronic version. 



8 



initial values 
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Figure 1: Evolution of the 
BI universe for a positive a. 




-1 -0,5 0.5 

ATAN{I) 



Figure 2: Evolution of the 
BI universe for a negative a. 
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0,9 



1,4 



Figure 5: 3D view of arctan(p), arctan(//) and arctan(T) for a > 0. 



The system of four time-depended variables T, //, p , T possesses integral curves 
in 4D space. The graphs of p{t) and T{t) are the projections from the extended 
space of dimension 4+1 on the plane with dimension 1+1. Therefore, the projec- 
tions of integral curves can intersect. In the figures the integral curves correspond 
to different initial conditions with similar other parameters. In the figures [5] and [6] 
the dot lines leads to the past, the solid ones to the future, while the point is the 
one where the initial condition of integration was given. 

There is also the possibility to assume that S is constant and the analysis de- 
velops in the aforesaid manner, practically interchanging a2 with . 

3.2.2 Assumption 2: R and S differ by a constant 

Let us now assume the case when R — S = r2 with r2 a constant. This assumption 
leads to 

a2 = C3a3T'"2/^ (3.31) 

which together with (12.1 II ) gives 

T = To + (3J/r2)f. (3.32) 

Thus we find T to be a linear function of t. It should be noted that analogical 
results occurs when the BI Universe is filled with stiff fluid. On the other hand 
from (12.91) one finds 
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-0.1 04 
ATANCH) 



0,9 



1,4 



Figure 6: 3D view of arctan(p), arctan(//) and arctan(T) for a < 0. 



ai = (l/c3)a^^Ti"''2/3. (3.33) 
In view of (l332l) . (l333l) and (I3ll1) from (I2l2l) one finds 

P = - - , n 4 ■ (3.34) 

Finally, inserting (l334l) . (l333l) into (I2l3]) we find 

a3 = c4T(«(3-2'-2)-3+.2)/6(2a-l) ^ (335) 

with C4 another constant. 

Thus the Einstein system of equations has been completely solved . 

3.2.3 Assumption 3: R and S are proportional 

Let us now assume the case when R = r^S with being some constant. This leads 
to the relation 



= ci3, (3.36) 



a2 
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with ci3 being some arbitrary constant. This relation together with (12.1 II) and 
(|Z9l ) gives 



ai = 4-+i)/3(i--)c^f(^--)TV3e^^^-, (3.37a) 



a2 = c^2^^)/^(^-^^^)c7//^(i-'-^)t1/3.^^^-, (3.37b) 



a, = 4--2)/3(i-.3)^-V3(i-.3),i/3^f^/f_ (337^) 

The system analogue to (I3.30|) is: 

t = 3//T, (3.38a) 
T = -, (3.38b) 

V3al-r3T a / 

This system is quite similar to (13.301 ) and choosing suitable parameters we find 
solutions which look as the ones illustrated in the previous graphs. 



4 Summary and outlook 

We have offered an investigation of an anisotropic cosmological BI model. Hav- 
ing in mind the complexity of the model we used some tractable assumptions 
regarding the parameters entering the model. For different assumptions regard- 
ing the shear anisotropy parameters (I3.25al) and (I3.25bl) we get interesting models 
deserving the study. The analytical results are supplemented with numerical and 
qualitative analysis describing the evolution of a BI Universe for different values 
of the parameters. 

In our further studies fTS*] we should like to see how the model isotropises 
at late times. Any realistic model must lead to isotropisation necessary for com- 
patibility with standard cosmological models at late times and in agreement with 
current observations. Also it is important to know how stable is the model to 
perturbations of the parameters. 
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